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We present a method for calculating interatomic exchange interactions within density-functional theory (DFT) 
in the full-potential linearized augmented plane wave (FLAPW) scheme. Our approach is based on total-energy 
calculations of spin-spiral states in reciprocal space employing the generalized Bloch theorem and a Fourier 
transform in order to obtain the real-space interactions, including the case of more than one atom in the unit 
cell. We discuss applicable symmetry relations that allow for a reduction of the computational time. We examine 
under which conditions the force theorem holds in close-packed and open systems, in particular regarding the 
treatment of the magnetization in the interstitial region. By including the exchange interactions into a classical 
Heisenberg Hamiltonian, we study the thermodynamics at finite temperatures within a Monte Carlo method. We 
specifically discuss necessary corrections in the case that the magnetic moments of one of the atomic species 
are soft or induced by strong moments of the neighboring atoms. We argue that a quadratic dependence of 
the energy on the soft local moment is a good approximation in many such cases, and for this situation we 
derive exact results which allow an elimination of the soft moments from the thermodynamic calculations by 
considering a renormalization of the strong-moment interactions. We prove that these renormalized interactions 
are temperature-independent. 
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I. INTRODUCTION 

In recent years there has been increasing activity in the pre- 
diction of high-temperature magnetic properties of solids, es- 
pecially regarding critical magnetic transition temperatures. 
The theoretical approach is founded in many cases on two 
assumptions: (i) that the magnetic excitations of a system 
can be phenomenologically described within a classical or 
quantum Heisenberg model, and (ii) that the exchange pa- 
rameters entering the model (i.e., the excitation energies) can 
be microscopically derived from total energy results of e.g. 
density-functional theory calculations. This "magnetic multi- 
scale modelling" has proven succesful in many cases, includ- 
ing itinerant elemental ferromagnets such as Fe and Cor^ 
localized-moment systems such as Gd,'''- magnetic alloys as 
NiMnSb or Co2MnSi^ii or diluted magnetic semiconductors 
as (Ga,Mn)Asi^ In these systems, the estimated Curie tem- 
peratures Tq are within 10-15% of the experimental values, 
showing that the approach is reliable for practical purposes. 

The derivation of magnetic excitation energies from 
density-functional calculations requires additional assump- 
tions, since density-functional theory is, in principle, valid 
only for the description of the ground state of many-electron 
systems. Mainly, one relies on an adiabatic hypothesis, which 
conjectures that the slow motion of low-energy magnetic ex- 
citations can be decoupled from the fast motion of intersite 
electron hopping, so that the local electronic structure has time 
to relax under the constraint that a magnon traverses the sys- 
tem. Then the magnons are regarded as practically static or 
"frozen" objects, and constrained density functional theorjiiS 
is employed. This adiabatic approximation, together with the 
realization that the local moments persist until above the Curie 
temperature, constitutes a paradigm which has proven fruitful. 
Quite a few theories and methods have been based on it, aimed 
at the calculation of thermodynamic quantities, including Tc, 



by harvesting the model parameters from first principles cal- 
culations and working out the thermodynamics within Monte 
Carlo methods or other suitable approaches. 

The adiabatic approximation in spin dynamics was dis- 
cussed by Gyorffy and co-workers in the development of a 
mean-field theory of magnetic fluctuations." Later Antropov 
et alpi2, and Halilov et alji^ derived equations for adiabatic 
ab-initio spin dynamics. Further elaboration came from Niu 
and collaboratorsii and Gebauer and Baroni,'^ who showed 
that the Berry curvature enters the adiabatic dynamics equa- 
tions; they also demonstrated mathematically how the Born- 
Oppenheimer method can be generalized for adiabatic spin 
dynamics without the requirement that there should be a large 
mass difference defining the time-scale difference. 

Based on the adiabatic approximation, a proper 
parametrization of the excitation spectrum is called for. 
In many, but not all, ferromagnetic systems the magnitude 
of the moment magnitude is relatively robust under modest 
rotations, so that mainly the pair exchange constants entering 
a Heisenberg model are required. From the point of view 
of methodology, two approaches are commonly used for 
the determination of the pair exchange constants. The first 
is based on multiple-scattering theory and Green function 
methods, frequently in the approximation of infinitesimal 
rotationSfi^ It is widely used used in methods where the 
Green function is available j ' --i^i ' 

The second approach, which we follow in this work, is 
based on reciprocal-space calculation of excitation energies 
by constraining the system to spin spirals of given wavevec- 
tors q. This requires non-collinear calculations, restricted 
to the primitive chemical unit cell (in the absence of spin- 
orbit coupling) by virtue of a generalized Bloch theorem;^ A 
subsequent integration in q-space (in form of a back-Fourier 
transformation) yields the pair exchange constants. Although 
the force theorerr^ is in principle not necessary here, prac- 
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tically it is often used in order to avoid a the numerical load 
of a self-consistent calculation for each vector q. This sec- 
ond approach is well-suited for electronic structure methods 
which are based on Hamiltonian diagonalization rather than 
Green functions, and has been developed, for example, for the 
augmented spherical wave (ASW) method^^^ or the LMTO 
method. ' - The results of the two approaches to the calculation 
of exchange interactions agree well, if the density-functional 
calculations are done within the same electronic structure 
method.-" 

This paper contains two major parts. In the first part. Sec- 
tions and |III1 we present a formalism for the calculation of 
pair exchange constants within the full-potential augmented 
linearized plane wave (FLAPW) method.— The approach of 
spin spirals and inverse Fourier transformations is used to ar- 
rive at formulae for the pair exchange constants in the case of 
one or more magnetic atoms per unit cell. Accuracy tests are 
shown concerning the use of the force theorem at finite rota- 
tions; we also address the problem of subtracting the contribu- 
tion of the magnetization in the interstitial. As a test, we calcu- 
late the Curie temperature of certain compounds by a Monte 
Carlo method. In the second part. Section HVl we discuss a 
way to parametrize the energetic contribution of longitudi- 
nal changes in the atomic spin moment so as to include them 
in a Monte Carlo method in the case that we are faced with 
a compound containing strong-moment and induced-moment 
sublattices. We obtain a scheme for the study of temperature- 
dependent magnetic properties and derive equations that allow 
a simplification of the computational method and a reduction 
of the computational cost under certain frequently-met phys- 
ical conditions; in particular we show how the weak-moment 
degrees of freedom can be exactly eliminated in favour of 
renormalized, temperature-independent strong-moment inter- 
actions. An application on the half-Heusler alloy NiMnSb is 
presented. 



II. AB-INITIO CALCULATION OF HEISENBERG 
EXCHANGE PARAMETERS 

Adopting the adiabatic approximation, the magnetic inter- 
actions are modelled by a classical Heisenberg Hamiltonian. 
The part of the total energy due to these interactions is then 
obtained from the expression 



E = 
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where R,na(np) = Rm{n) + T„(/3). Here, -R,„(„) are the 
lattice vectors and T^i^p) are the basis vectors specifying the 
positions of the atoms within the unit cell. Mmaini3) are 
the atomic magnetic moments at the sites Rma{nf!)^ while 
J{Rr,ia, Rnp) IS the cxchangc coupling constant for the pair 
of atoms situated at these sites, and is the quantity to be cal- 
culated. The summations over the indices n, m are carried out 
over all lattice vectors, and the ones using indices a, (3, over 
all the atoms in the unit cell. The factor 1/2 takes care of the 
double counting and the on-site term (Rma — Rnp) is left 
out. 



The constants J {Rma, Rnp) contain the information about 
the inter-site interaction due to the exchange coupling. The 
knowledge of these exchange interactions is essential for the 
description of thermal excitations in magnetic solids and their 
deriving from ab-initio calculations is the core problem in the 
attempt to describe the system with the Heisenberg Hamilto- 
nian. The correspondence between the ab-initio theory and 
the Heisenberg model is estabUshed by using the ansatz 
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which follows from Eq. ([Hi, as a defining relation of 
J {Rma, Rnp) within density-functional theory. Here, SMna 
and SMn0 are to be understood as small differences with 
respect to the direction only, not the magnitude. I.e., an 
appropriate energy functional (usually within the local den- 
sity or generalized gradient approximation) i?DFT [p, "m] of 
charge density p{r) and magnetization density m(r) is used 
in the place of E in Eq. (|2]i. When evaluating (|2]i, it is as- 
sumed that the exchange-correlation field is collinear within 
each atom, so that the derivative with respect to atom-cell 
integrated moments Mma = Jcdima ''^{f) d^''' is meaning- 
ful. The intra-atomic collinearity is an approximation, jus- 
tified by the energetic dominance of the moment formation 
(usually of the order of eV) compared to the formation of 
ferromagnetic order (usually of the order of less than 0.1 
eV). From these comments it is also evident that we do not 
require that the local moments are quantized either in the 
form M = y^S{S + 1) or = S" with S inte- 
ger or half-integer, as, e.g., would be the case in ferromag- 
netic semiconductors.-" Rather, Eq. ^ represents the lowest 
term in an expansion of the total energy in terms of the mag- 
netization direction, neglecting longitudinal enhancement or 
suppression of the moments, and Eq. ^ represents a classi- 
cal Heisenberg model, valid after local quantum spin fluctua- 
tions have been averaged out (see, e.g., the discussion in jlsll ). 
There are known cases when these approximations are insuf- 
ficient, and we discuss such a case in Section llV] 

Next, the collective transverse magnetic excitations are ap- 
proximated by static spin spirals, the energy of which is calcu- 
lated within the non-collinear FLAPW method.-'' In the case 
of a spin spiral with wave vector q, the azimuthal angle of the 
magnetic moment of an atom at the position Rna is given by 
fna = Q ■ Rna- The magnetic moment of an atom at the 
position Rna is 



Mna = Ma 
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sin 9a sin(q • R 

cos 9a 




(3) 



where 9a is the so-called cone angle, a relative angle between 
the final and initial direction of the local magnetic moment 
(here chosen along the z axis; this choice does not limit the 
generality in absence of spin-orbit coupling), and (pa an even- 
tual phase factor, also called phase angle. Taking advantage 
of the translational invariance we define R = ii„ — Rm and 
Tap = Ta — Tp, whcncc Eq. ^ becomes 
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Here, the energy E isa function of the spin-spiral vector q, as well as of the cone and phase angles of the magnetic moments on 
all the atoms of the unit cell. The dependence on these angles is collectively expressed by Q for the set of all cone angles {9a} 
and by $ for the set of all phase angles {i/iq,} in the argument of E. To account for the condition Ta ^ Tp — R under which the 
sum in Eq. (jUi is conducted (and is indicated by a prime), from now on we set J(tq,, Tq) = 0, for all the atoms a in the unit cell. 

With the aim to obtain the exchange interaction constants J(Ta , Tp — R) at the minimum of computational expense, we define 
in the following a set of expressions which are evaluated computationally. We first define the Fourier transform 
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It is straightforward to show that with the use of this Fourier transform, Eq. (HI becomes 
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A. Symmetry Relations 



B. Calculational Scheme 



Starting from the condition that J{Rma, Rnp) are real and 
symmetric and from the definition of the Fourier transform 
Jaf}{q) (Eq. 12]), several useful symmetry relations of Jaf}{q) 
can be derived (valid for each q vector): 

1. Jap{q) = Ji3a{-q) 

2. Tie [Jap{q)]=ne[Jap{-q)\ 

3. Im [Jai){q)] = -1m [Jap{~q)] 

(3a) Im[Jap{0)] =0 
mim[Jaa{.q)]^Q 

4. Jap{Cq) = Ja'13' (q), where (7 is a crystal point group 
symmetry element and (a' , f3') are the equivalent sites 
in the unit cell to {a, f3) via the action of the symmetry 
element C*"^, i.e., C~^Ta ~ Ta' + R for some lattice 
vector R (and analogously for /3). 

A very important consequence of Symmetry Relation 4 is that 
the q vectors can be sampled from the irreducible wedge of 
the Brillouin zone while Jap{q) in the rest of the Brillouin 
zone can be obtained by the symmetry transformations. In 
case that the crystal possesses the inversion symmetry / and 
if Tq, — iTa = 2Ta and Tp — iTp = 2Tp are both lattice 
vectors, then from Symmetry Relations 3 and 4 it follows that 
all Jap{q) are real. Moreover, due to Symmetry Relation 1, 
even if the system does not possess the inversion symmetry it 
is not necessary to make two separate calculations for q and 

-q- 



To develop a scheme for the calculation of the Fourier trans- 
forms Jap{q), we distinguish two different cases in the calcu- 
lational setup. 

•Case 1: All the atoms in the unit cell are ordered in the 
collinear ground state, except for atom /i. Its magnetic mo- 
ment is tilted by the cone angle 6*^ and the spin spiral running 
through the system will affect only the magnetic moments sit- 
uated on the atoms of the same kind as /i. In short: 

0^ + 

61a = 0, VAt^^ 



With the use of the symmetry relations for the coefficients 
Jap{q), from the total energy expression ^ one obtains 



A/2 sin^ 



(7) 



•Case 2: 

This case will appear only if there are two or more mag- 
netic atoms in the unit cell. Keeping the rest of the magnetic 
moments parallel, the magnetic moments on atoms /i and i/ 
are tilted by cone angles and di,, respectively, so the spin 
spiral running through the system changes the orientation of 
magnetic moments on both of these atoms: 
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As we have seen, if the system does not possess inversion 
symmetry, the coefficients J^^(q) are complex for ji ^ v. 
Their real and imaginary part can be obtained as 
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where E^''{q, 0) and E'^'^{q, ^) denote the total energies in the presence of a spin spiral defined with the wave vector q and the 
difference of the phase factors 0^ — = and ^ respectively. 



C. Brillouin Zone Integration 

We have established that the Fourier transforms J^^ (q) can be obtained from the differences in total energy between the 
states having specified magnetic configurations. Armed with Eqs. (|7][8]) and (|9]) we are now ready to calculate the Heisenberg 
exchange coupling constants, J(t^, t^ — R). First, however, one has to take into account that from the Eq. ^ it is only possible 
to calculate the difference J^^{q) — J^^ (O), but not the coefficient J^^{q) alone. This problem can be easily circumvented by 
introducing the coefficients J^u{q), defined as 



(10) 



Also, for simplicity, the non-zero cone angles can in all calculations be taken to have the same value 0. Eqs. (|2l[8]l and Q can 
now be re-written as 



7^e 



Et^iq) - £;^(o) 

ii/,2 sin^ e 



M^A/„ sin^ e 2JU' 
M^Af^ sin^ e 



2M.r 



Ue 



(11) 



The final step will be a simple back-Fourier transform. Using 
Eqs. (|5]l and ( fTOl i it is easy to see that 

yBZ jvbz 

Finally, we note that from the definition ( fTOl i it is clear that 
Jpy(q) satisfies the same symmetry relations as the coeffi- 
cients J^^{q). 

The described calculations can be time consuming, since 
they involve the determination of small energy differences 
(typically of the order of a few mRy). Due to the oscilla- 
tory phase in Eq. (fT2l) . the appropriate size of the q point 
set increases with the distance between the two atoms for 
which the interaction constant is being calculated, the grid 
fineness being basically determined by the inverse of the quan- 
tity \t^^ — R\ that enters the exponential in Eq. In the 
spirit of the Nyquist-Shannon sampling theoremj^ and for 
\R\ ^ |T7ti/|, the g-grid fine spacing in the direction of R 



should be at most half the value 2tt/\R\ . Additionally, suf- 
ficient accuracy requires larger plane-wave basis sets and a 
finer fc-point grid compared to a simple ground-state calcula- 
tion. A rule of thumb for increased accuracy is that, given a 
q-grid the fe-grid should be twice as fine per spatial dimension 
of the lattice in order to avoid spurious oscillatory behaviour 
of period of the grid-spacing Sk in J^.i,{q)- A self-consistent 
calculation of all energies needed here is computationally very 
demanding. Fortunately, in many cases the spin spiral can be 
considered a small enough perturbation that Andersen's force 
theorem2i can be used to calculate the energy differences. We 
discuss this in the following subsection. 



D. Test of the Applicability of the Force Theorem 

The magnon energy is calculated as the difference between 
the total energy of the system with a spin spiral (this is an ex- 
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cited state), and the ground state, which is ferromagnetic in 
the systems under study here. A self-consistent calculation of 
the spin spiral requires use of a constraint, in the form of an 
external spiraling magnetic field, which forces the magneti- 
zation to take the form given in Eq. (O. On the other hand, 
an application of the force theorem requires a position depen- 
dent rotation of the exchange-correlation field so that its 
direction acquires the form i.e., 

(sin 6a COs(q • Rna + 0q) \ 
sin 0„siii(q-ii„„ + , (13) 

cos 6a I 

where Ba'^ = {V^^ — V^^) is the self-consistent exchange- 
correlation field of the collinear calculation at atom type a 
(V^a is the exchange-correlation potential dependent on spin 
a =t, i). In the FLAPW method, this rotation is applied in the 
muffin-tin spheres, i.e., non-overlapping spheres around the 
atomic nuclei where the potential is expanded in radial func- 
tions and spherical harmonics, to be contrasted with the inter- 
stitial space between these spheres. Using the field of Eq. (fTsl l 
in the Kohn-Sham equations yields a (non-self-consistent) 
sum of eigenvalues of the occupied levels; according to the 
force theorem, the magnon energy is approximated by the dif- 
ference of this sum to the sum of eigenvalues of occupied lev- 
els in the self-consistent ferromagnetic ground state. 




(0,0,q) [ic/a] 



FIG. 1: (color online) Comparison of the force theorem with self- 
consistent calculations for NiMnSb. Left: a dispersion curve of a 
spin spiral along the [001] direction for a cone angle 6 = 30°. Right: 
spin spiral energy vs. sin^ 6 for a short-wavelength spiral along the 
[001] direction. 



to be a reasonable choice, since the energy differences are not 
too large for the magnon to stop being a small perturbation, 
but are also not too small so that one would have to employ 
a very large basis, or k-point set. This cone angle was used 
in the calculation of the Heisenberg interaction constants J 
presented in Section Hill 

We close this section with a comment on a spurious effect 
of the exchange-correlation field in the interstitial region when 
the force theorem is used. In a force-theorem application, 
the trial exchange-correlation potential in the interstitial is a 
smooth periodic function. However, in a self-consistent spin- 
spiral calculation, the resulting exchange-correlation potential 
is in many cases a much less smooth, although still periodic, 
function.-'. Then the interstitial magnetization in the force- 
theorem calculation can cause a serious overestimation of the 
spin-spiral energy. This spurious energy contribution can be 
circumvented by setting the magnetic part of the exchange- 
correlation potential in the interstitial to zero. Depending on 
the volume filling of the touching "muffin-tin" atomic spheres 
of the system, the spurious energy can be considerable, be- 
coming larger for open systems. In Fig. we show a calcu- 
lated example for NiMnSb. The self-consistent spiral energy 
and the force-theorem energy (here normalized to sin^ 6) are 
practically identical, if we set = in the interstitial in 
the force theorem calculation (in the self-consistent calcula- 
tion, B^'^ is non-zero also in the interstitial); if this correction 
is not used, then the spiral energy is strongly overestimated. 
Figs. Ilb-d show the same for Fe2MnSi in the full-Heusler 
structure, FeMnSi in the half-Heusler structure, and FeSi in 
the zincblende structure. All three structures are based on an 
fee lattice with a lattice parameter aiat = 5.663 A, but differ 
in the number of atoms in the unit cell and therefore in the 
volume of the interstitial region. Fe2MnSi has the smallest 
interstitial region, and evidently the interstitial magnetization 
has almost no effect; both force theorem calculations prac- 
tically coincide with the self-consistent result. For FeMnSi, 
however, the interstitial volume is larger and the spiral energy 
is strongly overestimated in the force theorem calculation, if 
the interstital magnetization is not set to zero. The effect is 
even stronger for FeSi, where the interstitial volume is largest. 



The approximation is expected to be better for smaller per- 
turbations, i.e., for smaller cone angles and/or smaller magnon 
wavevectors |q|. As a test, in the left panel of Figure [T| the 
dispersion curve is shown for a spin spiral in NiMnSb, de- 
fined by a cone angle 6 = 30°, and a spin-wave vector along 
[001] direction. The force-theorem and self-consistent cal- 
culations agree rather well. The right panel of Fig. [T] shows 
the dependence of the magnon energy on the squared si- 
nus of the cone angle (sin^ 6) for a fixed spin-spiral vector 
q = (0,0, l)7r/aiat. We see a a maximal deviation of the 
order of 7-8% for the unfavourable case of 6 ^ 90° and 
q ~ 7r/aiat, while the deviation starts becoming visible at 
cone angles larger than 0„ia.T ~ 50°. 

A general conclusion is that if one wants to use the force 
theorem to obtain the spin-spiral dispersion within the whole 
first Brillouin zone of the crystal, a cone angle 6 — 30° seems 



III. EXCHANGE INTERACTION PARAMETERS AND 
CURIE TEMPERATURE 

Following the prescription of Sec. HI] we calculated Heisen- 
berg exchange interaction parameters (J^; i,j are abbrevia- 
tions of ma, n(3) of Co2MnSi and NiMnSb, shown in Fig. [3] 
and used a Monte Carlo method to obtain an estimate of the 
Curie temperatures for these compounds. The calculations for 
both compounds were performed within the generalized gra- 
dient approximation on a 4096 fc-point mesh with 2744 q- 
points in the full Brillouin zone. The planewave cutoff was 
^max = 3.8 au^^. The convergence was checked with respect 
to the above parameters. 
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FIG. 2: (color online) (a): Spiral energy E{q, 9) as function of sin^ 6 in NiMnSb for the cases of (i) a self-consistent calculation (including 
the self-consistent, spiraling polarization of the interstitial) (ii) a force-theorem calculation including the polarization of the interstitial, and 
(iii) a force-theorem calculation excluding the polarization of the interstitial. Evidently, case (ii) stongly overestimates the spiral energy, since 
the interstitial polarization acts as an effective magnetic field. Case (iii) on the other hand is a good approximation to the exact result (i). The 
spiral wavevector is q = (0, 0, 0.15) (vr/aiat). (b,c,d): Same as in (a) but for Fe2MnSi, FeMnSi, and FeSi, in fuU-Heusler, half-Heusler, and 
zinc-blende structures, respectively, to demonstrate the effect of increasing the volume of the interstitial. The structures correspond to an fee 
lattice with four, three, and two atoms per unit cell, respectively. All three were calculated in the same lattice parameter (aiat ~ 5.663 A) and 
at a spiral vector of q = (0, 0, 0.25) (vr/aiat). 
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In Co2MnSi, both Co and Mn atoms have strong and sta- 
ble magnetic moments, whose interaction is described with J 
parameters depicted in Fig.[3h. In NiMnSb, though the mag- 
netic moment of Ni is small and actually induced by the Mn 
surrounding atoms, both Ni and Mn were treated as magnetic 
atoms and the parameters of Mn-Mn, Ni-Ni, and Mn-Ni in- 
teraction were calculated (Fig. [Sti). As will be discussed in 
Sec. IIVI this treatment gives an insight into the thermal be- 
haviour of the Ni sublattice and with the use of an extended 
Heisenberg model more useful information can be obtained. 

Co2MnSi and NiMnSb are half-metallic ferromagnets, i.e., 
the density of states (DOS) in one spin direction (here ma- 
jority spin) is metallic, while in the other spin direction the 
DOS shows a band gap around the Fermi level. For a non- 
half-metallic ferromagnet, the interaction constants follow 
the characteristic Ruderman-Kittel-Kasuya-Yoshida (RKKY) 
asymptotic behaviori 
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where kp is a Fermi wave vector direction such that the as- 
sociated group velocity is parallel to Ri — Rj and ip de- 
notes a phase factor The exchange interaction according to 
Eq. (fTSl l has an oscillatory character with an envelope decay- 
ing as \Ri — Rj\^. On the other hand, in a half-metal there 
are no states at the Fermi level in the minority band, thus the 
relevant Fermi wave vector is imaginary, kp = iKp, cor- 
responding to decaying states. Hence, one obtains an expo- 
nentially damped RKKY behavioral In both cases shown in 
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FIG. 3: (color online) Heisenberg exchange interaction parameters of 
Co2MnSi (a) and NiMnSb (b) as a function of the distance R between 
the atoms (in units of the lattice constant aiat). The lines are guides 
to the eye. 
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FIG. 4: Heisenberg exchange interaction parameters of fee Co as a 
function of tfie distance R between the atoms (in units of the lattice 
constant a). The line is a guide to the eye. 
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by the peak in the temperature-dependent, static susceptibil- 
ity. As the simulation supercells are rather large, finite-size 
corrections to Tq are small. 

We also show Monte Carlo magnetization curves for 
Co2MnSi in Fig. |5] The classical Heisenberg Hamiltonian, 
Eq. ([T]i was used to model the systems, assuming that the 
magnetic moments can change their orientation, but not their 
length. NiMnSb is discussed in more detail in Sec. |IV] In the 
calculations the Metropolis method was employed. A super- 
cell of 2048 unit cells (4096 magnetic atoms) was used for 
NiMnSb and one of 2197 unit cells (6591 magnetic atoms) 
for Co2MnSi; interactions to neighbors up to a distance of 
four lattice constants were included. For each temperature the 
number of Monte Carlo sampling events was 5000, after al- 
lowing an initial relaxation time of 1000 steps and taking one 
sampling event every 10 sweeps of the lattice. 

The magnetization curves (Fig.|5]) do not go to a sharp zero 
at Tc, but rather have a tail, as a result of the finite supercell. 
The peak of the susceptibility is, on the other hand, rather 
sharp and its position can be used to determine the Curie tem- 
perature. From the positions of these peaks, we estimated for 
CosMnSi (Fig.m left) rc=980 K (experimental^ 985 K), and 
for NiMnSb (see Sec. |IVl) Tc=870 K (experimental^? value 
being 730 K). Note that here, for NiMnSb, the Ni moment 
was also taken into account within the Heisenberg model (but 
see also the discussion in Sec. lIVI in particular llV Cb . For fee 
Co we find a Curie temperature of 1200 K, while the experi- 
mental result is 1403 K (the high-temperature stable phase of 
Co is fee). 



IV. CONSIDERATION OF LONGITUDINAL MOMENT 
FLUCTUATIONS 



FIG. 5: (color online) Magnetic moment AI{T) per unit cell (blue 
full line), partial moment of the Mn and Co sublattices (green dashed 
lines), and susceptibility x(T') (red full line) as functions of the tem- 
perature T, for Co2MnSi. The peak of the susceptibility at 980 K 
indicates the Curie temperature (experimentally found to be 985 K). 
Note that the susceptibilities are shown in an arbritary scale. 



Fig. [5] we notice this very fast decay of the interaction param- 
eters with the distance between the atoms. For a comparison, 
in Fig. 151 the exchange interaction parameters of fee Co are 
shown. Obviously, the decay here is much slower. 

Once the exchange interactions are known, the effective 
Heisenberg model can be solved for the Curie temperature. 
While the mean field approximation is computationally the 
simplest method to this task, it is known (and has been shown 
in practice) that the resulting Tc is overestimated. The ran- 
dom phase approximation (also known as Tyablikov approxi- 
mation), on the other hand, is rather accurate."*^ The most ac- 
curate, but numerically more expensive, way to calculate Tc 
of a classical Heisenberg model is the Monte Carlo method, in 
particular when taking advantage of the cumulant expansion 
to account for the finite size of the simulation supercells. In 
this work we applied the Monte Carlo method, locating Tc 



So far we have discussed a formalism and examples of cal- 
culation of exchange parameters for a Heisenberg model, as- 
suming that the magnitude of the local magnetic moments is 
rigid. However, as is long known, certain systems are weakly 
magnetic and the moment formation takes place at a rela- 
tively low energy scale, comparable to the magnon energies. 
Moriya's quantum mechanical spin fluctuation theory^ al- 
ready factors in these effects. An introduction of a classical 
Hamiltonian including longitudinal and transverse degrees of 
freedom has been done by Uhl and Kiibler?- by parametrizing 
ab-initio total energy results; thermodynamic quantities were 
then obtained within a mean-field approach which coupled 
longitudinal and transverse degrees of freedom^SMS More 
elaborate approaches by a classical fit of the energy to density- 
functional results, including transversal and longitudinal de- 
grees of freedom, together with Monte-Carlo or classical spin- 
dynamics calculations, were presented e.g. by Ruban et al.r^ 
Ma and Dudarev,?* or Derleti^ 

Particularly interesting are compounds of a strongly mag- 
netic and a weakly magnetic subsystem, where the weak mo- 
ments are cannot be treated as having rigid-magnitude, but 
on the other hand are large enough that they cannot be ne- 
glected. Examples are FePt, FePd, or NiMnSb, where the 
strongly magnetic atoms are Fe and Mn, while the weakly 
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magnetic atoms are Pt, Pd and Ni. Then the assumption of 
rigid Heisenberg spins does no longer seem plausible. One 
way to circumvent this is to use only the strong magnetic mo- 
ments as independent variables, but with their pair exchange 
parameters renormalized by the exchange among the weakly 
magnetic atoms and by the enhanced susceptibility. Such an 
approach was, for example, used by Mryasov et alM^ to 
model the temperature dependence of the magnetic anisotropy 
in FePt. Sandratskii et a/4S also find improvement in the theo- 
retical results if the weak moments are not treated as indepen- 
dent, rigid Heisenberg spins, but rather fully adjust to the mag- 
netization of the strong-moment sublattice. There is, however, 
the point of view of treating the weak moments not as fully 
enslaved to their strong-moment environment, but as indepen- 
dent degrees of freedom whose fluctuation can affect the ther- 
modynamics. This point of view has been pushed for example 
in order to understand the antiferromagnetic-to-ferromagnetic 
transition of FeRh^^^ 

In this section we develop a method of treating such sys- 
tems, emphasizing the choise of the energetically relevant 
magnetic degrees of freedom at temperatures up to the crit- 
ical temperature, as well as the practical implementation of 
these degrees of freedom to models that can be solved via 
Monte Carlo simulations. The additional challenge, compared 
to the theory developed in the previous sections, stems from 
two facts. First, a treatment of the weak moments as indepen- 
dent degrees of freedom must include the energy scale of the 
moments magnitude as well as their direction. Second, if the 
interactions of the strong moments are calculated as described 
earlier, by considering non-collinear magnetic configurations, 
then the calculated energies include a spurious contribution 
because the magnitude of the weak moments changes during 
the application of the non-collinear constraint. This spurious 
contribution must be accounted for However, at the end we 
see that, in certain commonly occuring situations, the degrees 
of freedom of the weak-moment atoms can be eliminated 
in the statistical-mechanical calculations by using renormal- 
ized exchange parameters, while the statistical averages of the 
weak moments follow directly from the averages of the strong 
moments. We show that this is an exact result if the weak 
degrees of freedom give only quadratic contributions to the 
energy. Here we make a step-by-step approach for the special 
case where the weak moments do not interact with each other, 
while we present the more general case when they interact in 
the Appendix. 



A. General approach 



The magnetically strong atoms are the main carriers of 
spin moment that also acts as an effective field po- 
larizing the magnetically weak atoms. The polarization of 
the latter is strengthened by an enhanced local susceptibility, 
which (neglecting complications from the band-structure) is 
X = XPauii/(l - /XPauli) (with xpauli the Pauli paramag- 
netic susceptibility and / the exchange integral). Let us denote 
by M the rigid-magnitude local moment of the magnetically 
strong atoms and by fi the supple local moment of the mag- 



netically weak atoms. A simple expression that parametrizes 
the energy in terms of the local moment fi and the polarizing 
field H"^ is 



(16) 



It is imphed that J?°« = «;E„e„cighb + H'^""'^ where 
Mn is the summed moment of the neighboring atoms induc- 
ing a polarization, k is a phenomenological parameter which 
encapsulates all microscopic processes (in particular electron 
hoppings) that couple /j, to Mn, while is an external 

magnetic field that we henceforth we set to zero4^ In the ab- 
sence of a field H'^^ Eq. ( fTSI l contains only even powers of fi, 
because symmetry requires E\oc[0, fi] = £'ioc[0, — /x]. 

In Eq. (fTSI l we always must have & > 0. The case of spon- 
taneous polarization is described by a < and 6 > 0, while 
the case of induced magnetic moments is described by a > 0. 
In this case, the on-site susceptibility is 



X - d^i/dH'^^'lH..^, - l/(2a) 



(17) 



Thus the enhancement of the susceptibility is contained in the 
parameter a. 

Eq. ( fT6b , if the fourth-order term is negligible, can be writ- 
ten as 



1 

2a' 



1 



-H'^tt if 6 = 0. 

/ 4a 

(18) 

To elaborate on these ideas we use NiMnSb as a concrete 
example. In this case the Mn subsystem is magnetically strong 
and the Ni subsystem magnetically weak. We have also found 
that in this case it is a good approximation to set 6 = (as 
we show below), which we adopt. Combining Eq. (fT6b at 6 = 
with the Heisenberg-model energy expression for the Mn 
subsystem results in the extended model Hamiltonian 

/GNi n{l) 

where ^.,^(1) '■^^ °f moments of the Mn near- 

est neighbors of the l-lh Ni atom, n{l), and Mi and /x; refer 
to the Mn and Ni moment respectively. The superscript "6" 
in the Mn-Mn exchange interaction is placed in anticipation 
that it does not coincide with the quantity calculated within 
the method presented in Sec. m but with a "bare" quantity, as 
will be discussed below. It has been also assumed that the Ni- 
Ni interaction can be neglected, as there are no Ni-Ni nearest 
neighbor pairs, while from the calculations it becomes evi- 
dent that more distant Ni-Ni interactions are negligible (the 
case where the interactions among the weak-sublattice mo- 
ments are non-negligible is discussed in the Appendix). The 
above equation can be rewritten in a somewhat more conve- 
nient form, if we define (with an obvious notation) 



Mn-Ni 



K, nearest neighbours Mn-Ni 
0, otherwise 



•^Ni-Ni = 



(20) 
(21) 
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Then Eq. ( fT9] ) takes the form 

H = - 



J2 4m,-m,+ j2 

ijeMn zeMn,ieNi 



J2 

jGMn.iGNi iGNi 



(22) 



This form is appealing, as it contains a Heisenberg-Hke ex- 
pression among all magnetic atoms plus a local correction 
term that accounts for the longitudinal degree of freedom of 
the Ni moments. However, it is important to remember that 
the Ni moments included in the scalar products of the RHS 
can change in size as well as angle, deviating from the tradi- 
tional Heisenberg model. 

We can derive the parameters k and a from DFT calcula- 
tions as follows. From Eq. ( fTSl ) it follows that the equilibrium 
value of the Ni moment at T = and for zero external field is 
linearly dependent to the polarizing neigboring Mn moments: 



(23) 



where TVc = 4 is the coordination number of a Ni atom with 
respect to Mn neighbors. From this expression one can deduce 
the ratio k/q once M and have been calculated in the fer- 
romagnetic state by an ab-initio method. On the other hand, 
a can be determined by introducing into the DFT calculation 
a constraint on the Ni moment, either in the form of a lon- 
gitudinal magnetic field constraining the magnitude or in the 
form of a transverse magnetic field constraining the angle 9 of 
fi with respect to the moment directions at the Mn neighbors, 
but allowing the magnitude to relax to /icq(6'). The former 
method results in a parabolic (for b = 0) moment-dependence 
of the total energy, as suggested by Eqs. (I16I18I ). The latter 
method results in the following dependence, as can be easily 
found by an energy minimization of Eq. ( fTSl ) at a given 9: 



Eioc{9) - £;ioc(0) 



Aicq(O)cos0 (24) 

a(/xeq(e) - Mcq(O))' 



= a^oq(O)^ sin^ 6*. 



(25) 



If the extended Heisenberg model provides a good approx- 
imation to the energetics of the system, then Eqs. ( fTSl l and 
( |25] | should both yield a good approximation to the calculated 
DFT energy difference, when a constraint is applied on the 
Ni moment, and the extracted parameter a for the two cases 
should be approximately the same. We find that this is the 
case in NiMnSb (see Fig This allows us to extract the 
value of a from the DFT calculation and from this the value 
of K (equivalent to JMn-Ni") ^^1- * l23T l. Thus the ingredients 
of formula (l22l) are accessible. 

There remains, however, a correction to be made, connected 
to a "renormalization" of the Mn-Mn exchange constants 
which have been calculated within the spin-spiral method de- 
scribed in Section [III To clarify the problem we remind the 
reader that, when calculating spin spirals, one normally con- 
strains the direction of the moments but not their magnitude. 
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FIG. 6: (color online) Top: Dependence of the total energy on the 
Ni moment in NiMnSb under the constraint of an external field act- 
ing on the Ni atoms. Circles: calculated data (within DFT). Dashed 
curve: parabolic fit to Eq. (TsJ with a = 53mRyd//iB. Bottom: De- 
pendence of the Ni moment and the total energy on the constraining 
angle 6 tilting the Ni moment away from the Mn moment in NiMnSb. 
Blue squares: /ioq(^) (scale displayed on the left-side ordinate). Red 
circles: EuFTid) — -Edft(O) (scale displayed on the right-side or- 
dinate). The broken lines correspond to fits to Eqs. l l24b and {25}. In 
the energy fit, the last part of Eq. ( 125b was used with the prefactor 
a = 55.3mRyd//iB fitted to the low-angle contributions (6 = 15° 
and 30°). 



For the Ni atoms, where the rigid-moment approximation is 
not vaUd, the magnitude /i is reduced by the spin-spiral for- 
mation due to the canting of the neighboring Mn moments in 
different directions. This effect provides an additional, indi- 
rect energy contribution to the spin spiral and to the Mn-Mn 
interaction, compared to the case that the Ni-moments magni- 
tude would have been kept constant. It is the spurious contri- 
bution that we mentioned in the introduction to Sec. |IVl Let 
us call this contribution JMn-Ni-Mn- The calculated Mn-Mn 
interaction consists thus of two parts: 



"^Mn— Mn — "^Mn-Mn 



Jun- 



■Ni-Mm 



(26) 



with Jmii-Mii '^he sought-after bare interaction entering 
Eqs. ( |20b and (l22l) . while JMn-Mn is a renormalized inter- 
action that is probed by the spin-spiral DFT calculation. For 
more distant Mn atoms, which have no common Ni neighbor, 

"^Mn-Mn COincidcS with JMn-Mn- 

For the derivation of an expression, e.g. in the case of 
NiMnSb, giving t/Mn— Ni— Mm considcr Nc{ — 4) Mn atoms 
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as nearest neighbors of a Ni atom. Then the local-energy ex- 
pression ( fTSl ) becomes 



Eloc = ^'^ X! ' ^ 

n=l 

= _A y !^M„ • M„, - — iV, Af 2 (27) 



Here, n and n' run through the Mn moments. The second step 
follows under the condition that, in the DFT calculation, the 
Ni moment relaxes to the particular equilibrium value that is 
dictated by the neighboring Mn-moment directions, i.e., /j, = 
^ X^n From Eq. ( |27] i we obtain the indirect interaction: 



Jun- 



loc 



Ni-Mn 



5M„5M„, 



K 

2a 



(28) 



At this point it is important to note that, whether one 
chooses to work with Jmh-Mii or Jmh-Mh' depends on the 
choice of the degrees of freedom. If only the Mn moments 
are chosen as degrees of freedom, then Jmh-Mh has to be 
used. If, however, the Ni moments are also chosen as inde- 
pendent degrees of freedom, then Jmii-Mh used to- 
gether with J^n-Ni ~ latter case, according to the 
prescription leading to Eqs. ( l20b and ( l28b , we arrive at the ex- 
tended Heisenberg Hamiltonian ( |22] ) with the appropriate bare 
parameters in the nearest-neighbor coupling. 



B. Analytical elimination of weak degrees of freedom at T > 

In the previous subsection we discussed the bare and renor- 
malized parameters of the model arising from total-energy 
calculations of the constrained ground state. In this section 
we examine the case of thermodynamic quantities at T > 0, 
where it is not a priori obvious that the same renormalization 
is still valid if the weak moments are allowed to fluctuate. We 
find that, in the absence of fourth order terms in E\oc [b ~ 
in Eq. (fT6ll1, also at T > the weak moments can be elim- 
inated in favour of the same renormalized parameters as the 
ones appearing at T = 0. Our conclusion is based on an exact 
analytical integration of the weak-moment part of the partition 
function in the case 6 = 0. 

First we observe that the energy functional ( flSl l under 
the action of has the same quadratic form as the one 



with H'^^ = 0, only with the minimum shifted to fi^q = 
H°^/{2a). This simplifies the integration of the partition 
function. To see this, we first transform the Hamiltonian ( fT9l ) 
in such a way that the renormalized interactions Jij appear 
explicitly. We start by rewriting Eq. iT% as: 



H{{M,;fii}) 




I n{l) 



(29) 



The indices i, j run over the Mn atoms, the index / over the Ni 
atoms, and n{l) over the Mn neighbours of the l-th Ni atom. 
The last term is now split in an interatomic contribution and 
an on-site contribution: 



I n(l) 



I n{l),n'{l) I 7i{l) 

i ^ gc(^,,)M..M,+^^gM„^ (30) 



2 ^ 2a 



I n(l) 



In the last step we have converted die sum n'(i)n^n' 
to a sum over i, j, by introducing a combinatorial factor c(i, j) 
that counts how many common Ni neighbors the i-th and j-ih 
Mn atoms have. From the structure of NiMnSb it follows that 
= 1, if hj are nearest-neighbors in the Mn fee sub- 
lattice (i.e. if the distance between i and j is aiat/\/2) and 
cihj) = otherwise. Actually we have recovered the quan- 
tity JMn-Ni-Mn = / (2a) of Eqs. dZTb and (|28] |. Finally, the 
last term of Eq. ( l30b is just a constant which can be omitted. 
Using the definition (|26l ) and Eq. (|28] |. we combine the first 
term of the RHS of Eq. (|29l) with the first term of the RHS of 
Eq. (l30b to obtain the renormailized . Then the Hamilto- 
nian to be used in the partition function takes the form: 



I \ n{l} 



(31) 



n{l) 
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where the Heisenberg Hamiltonian involving only the Mn sublattice with the renormalized interactions, Hr = ~ \ Tliij Jij^i ' 
Mj, has been introduced for convenience as part of the full Hamiltonian. Note that there is no renormalization for the Ni-Mn 
interactions: either they should be completely omitted within Hr that includes only the Mn moments, or they should be given 
by Eq. (l20b in the generalized model. The spin-spiral result for Jmh-ni does not enter in our theory. 
The partition function is: 

.3,, ^3„ mMf,M) 




dQi ■ ■ ■ dflN I d fii ■ ■ ■ d fiN exp <^ J~T J ^^^^ 

dn,---dnNexp'>-^} I dVi ■ • ■ d-'fiN exp <; - ° ^' ~ ^^"^'^ I (33) 

dni- --dnNcxp^--;^} { ^^^^^ ] (34) 

where J dfli ■ ■ ■ dil.N denotes an integration over the Mn-moment solid angles. The integration over the Ni-moments 
/ d'^fj.i ■ ■ ■ d^i-iN contains only the exponential of a complete square and has been analytically integrated to (nkBT/a)'^^^^, 
i.e., it corresponds to the partition function of 3N uncoupled harmonic degrees of freedom and is independent of the value of 
Mn- Zr is the partition function corresponding to the Hamiltonian Hr- The magnetization (Af) of the Mn sublattice also 
turns out to depend only on the renormalized partition function (and Hamiltonian). The same is true for any moment-moment 
correlation function of the Mn sublattice, {{MixY"'^' {MiyY'^'^y ■ - ■ (Mnz)"''^'), ct S (x, y, z), where the m„a > are integer 
exponents defining the order of the correlation function. To see this we consider explicitly 

((Mi,)"-(Mij,)"i-..(Mjv.)""=) = \ J dni---dnN{Mi..r'^ ■■■{MNzr"'exp{-nr/kBT} 

X / d>i • • • d>Ar exp j ^ ^ y (36) 



.^Jdn,--- dr!Ar(Mi.)"^== • • • (Mtv.)""^ cxpi-Hr /kBT} (37) 



where the Hamiltonian in the exponential has again been split 
in two terms according to Eq. ( [3Tl i and the integration over 
d'^/ii • • • d^HN has again been carried out analytically, can- 
celling out the factor (TrfesT/a)"^^/^ in the full partition func- 
tion. But expression (|37] | is just the expression for the corre- 
lation function of only the Mn sublattice with the renormal- 
ized Heisenberg interactions. Thus we see that, as far as the 
Mn-sublattice magnetization is concerned, one can use the ex- 
change interactions calculated by the spin-spiral method that 
are renormalized by construction, while neglecting the Ni mo- 
ments and the Mn-Ni interaction. If the Ni moments are to 
be included, then the bare interactions have to be used in an 
extended model, however, the Mn-sublattice properties in the 
two cases will be the same. Note, finally, that this exact re- 
sult is based on the fact that the exponent in the integration 
over d'^ni ■ ■ ■ d^^f^ contains a complete square, i.e., that the 
"harmonic" approximation, 6 = 0, is valid; if & 7^ and one 
proceeds to an elimination of the weak degrees of freedom, 
then the resulting renormalized strong-moment Hamiltonian 
can have higher-order terms (biquadratic, four-spin, etc.) with 
temperature-dependent parameters. 

A semi-analytical result follows in an analogous way also 
for the average Ni moment and fluctuation amplitude per 



atom: 



{^A = ^^c{M) (38) 

1 ^ 

= LlkBT+i—N,) l([yM,]2) 
a 2 \2a J N^^^ ^ ' 

1=1 

= -^fcBT+(f iV,)'(A/^) (40) 



I.e., there is a "harmonic" part and a part induced by the 
fluctuation of the Mn sublattice. The former is indepen- 
dent of the number of atoms, while the latter, (M^) := 
-^i]^) /-^' increases proportionally to the number of 
magnetic atoms in the system for T < Tq, as in normal 
Heisenberg systems4^ The longitudinal susceptibility xwi of 
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the Ni sublattice can be found in a completely analogous way: 



1 

2a 



(42) 



i.e., there is again a harmonic part and a part proportional to 
the Mn sublattice susceptibility. (In Eq. ( 1411 1. z is implied to 
be the direction of magnetization.) 



According to these results, as long as the approximation 
6 = holds, one can deduce the thermodynamic properties 
of the weak sublattice by merely a calculation on the strong 
sublattice, avoiding the extra numetical cost that a full Monte 
Carlo simulation entails. Actually this procedure can be car- 
ried out for higher-order correlation functions of the Ni mo- 
ments, {{fii^r'^-^ifiiyr^y ■ ■ ■ ifiNzr"') (a G (x,2/,z)), 
where once again the mia > are integer exponents defin- 
ing the order of the correlation function. In the resulting for- 
mula the integration over d^/i; can be carried out analytically, 
yielding a sum of correlation functions of the Mi of the form; 



J 



z 
1 



dQi--- dripfe ^ J fi>i • • • ■ ■ ■ iPNz)"^"' exp 



:r)<-"<^)"(£i:«» 

l^laexyz p^O ^ / ^ / nil) 

p even 



r 



rni{a)-p 



(43) 



where a change of variables + ^ ^na 

has been performed, the binomial expansion of {^la + 
^ J2n{i) -^na)™'" has been used and we have taken into ac- 
count that integrals of the type x^e""^ = y/^ip — 
l)!!/(2a)P/^ for even p and vanish for odd p [for the p = 
term we accept the convention (—1)!! = 1]. Due to the 
presence of products of AIna in Eq. ( |43] ) it is clear that this 
expression reduces to a sum of correlation functions of the 
Mi within the Hamiltonian Hr- The summations in this ex- 
pression are too involved to arrive at a general closed form, 
however, Eqs. ( l38l ). ( l40l i and (l42l) are special cases of appli- 
cation of this formula. Obviously, mixed correlation func- 
tions between strong and weak moments can also be reduced 
to strong-moment correlation functions of the renormalized 
Hamiltonian by first eliminating the weak moments following 
the same prescription. 

We should stress that choosing to eliminate the weak mo- 
ments in favour of renormalized interactions does not mean 
that the weak moments are physically less valid as degrees of 
freedom; to do so is merely a matter of mathematical or com- 
putational convenience, especially since efficient methods ex- 
ist for the calculation of thermodynamic quantities within the 
classical Heisenberg model. 



C. Calculations in NiMnSb 

We exemplify the above results with calculations in 
NiMnSb. First we establish which interactions can be ne- 
glected. To this end we performed calculations of the ex- 
change coupling parameters by the spin-spiral method. As it 
turns out, the Ni-Ni interactions are negligible. About the Mn- 
Mn interactions, Fig.[3}5 suggests that it should be enough to 



include up to 2nd neigbors, but this is misleading; we find by 
Monte Carlo calculations a 2nd neighbor approximation leads 
to an overestimation of approximately 100 K in Tc compared 
to the value including also more distant neighbors, therefore 
we include Mn-Mn interactions up to a distance of three lattice 
parameters. Among the Mn-Ni interactions (Fig.[3j3) only the 
nearest-neighbor coupling has some small influence, chang- 
ing Tc by a modest 20 K. Given the conclusions of the pre- 
vious subsection, the Mn-Ni interaction should be excluded 
for an estimation of Tq, since the Mn-Mn interactions are al- 
ready renormalized (in the corresponding spin-spiral calcu- 
lations the magnitude of the Ni moment was allowed to re- 
lax). The Mn-Ni interaction must be included if one wishes 
to extract information on the behavior of the Ni magnetiza- 
tion; then, however, the interaction type and strength has to be 
corrected, since the Ni atoms belong to a soft-magnetic sub- 
lattice. 

After deriving the exchange parameters by the spin- 
spiral method, utilizing the force theorem which according 
to Fig.|2^ is accurate enough for this purpose, we performed 
constrained-angle calculations for the Ni moments with full 
self-consistency. From the total-energy results shown in 
Fig. |6l together with Eq. ( 1251 ). we deduce a value of a = 
55.3 mRyd//x|. Given this, together with the ground-state 
magnetic moments of Mn M = 3.71 and of Ni /icq = 
0.26 /IB, Eqs. (|20l) and ^ yield a value of Jmh-ni = 
K = 1.94 mRyd//i|. This bare value Jmh-ni signifi- 
cantly different than the value of 0.92 mRyd//i| for the Mn- 
Ni interaction, which was derived from the spin-spiral cal- 
culation (Fig. [3j3). Finally, from Eqs. ( |26] | and ( l28b we can 
extract the bare value of the nearest-neighbor Mn-Mn inter- 
action: Jmii-Mii ~ 0.066 mRyd//i|, which shows a reduc- 
tion of about 1 /3 compared to the corresponding renormalized 
value. Note that the equilibrium moments have to be mul- 
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tiplied to these values if comparison is to be made with the 
energies shown in Fig.[3j3. Then one obtains Jun-Mn-^^^ — 
0.91 mRyd and JMn-Ni^Mcq = 1-87 mRyd; i.e., the Mn- 
Mn bare interaction is weaker than the Mn-Ni, which is not 
surprising, since the Ni polarization stems from a direct, 
nearest-neighbor exchange interaction with Mn. 




Temperature [K] 

FIG. 7: (color online) Top: Monte Carlo results on the magnetization 
of NiMnSb under various assumptions on the model and the interac- 
tions, (a-c): Only transverse fluctuations are taken, (a): Interactions 
as derived by force-theorem, spin-spiral calculations, including Mn- 
Mn and Mn-Ni. (b) Same as (a), but with bare Mn-Ni interactions 
■^Mn-Ni (Eq.[T9](. (c) Same as (b), but additionally with bare Mn-Mn 
interactions Jmh-Mii (given by Eqs. J26l28t ). (d-e): Also longitu- 
dinal fluctuations of the Ni-moment are allowed, (d) Bare Mn-Ni 
interactions, (e) Bare Mn-Ni and Mn-Mn interactions. Evidently, 
taking the Mn-Ni bare interaction alone leads to an overestimation 
of the excitation energies and Tc, which is corrected when the bare 
Mn-Mn interaction is also taken. Bottom: Magnetization curve and 
susceptibility of the Mn sublattice in NiMnSb calculated within dif- 
ferent models. Squares (magnetization) and circles (susceptibility): 
Only transverse fluctuations with renormalized Mn-Mn interactions, 
neglecting completely the Ni moments. Full line: Transverse and 
longitudinal fluctuations are allowed for the Ni atoms, only trans- 
verse for Mn, with bare Mn-Mn and Mn-Ni parameters; only the 
Mn sublattice magnetization and susceptibility is shown (i.e., cor- 
responding to the Mn part of curve (e) in the top panel). Note on 
the susceptibility units that l/ig/inR-yd = 4.254 x lO^^^ie/Tesla. 
The temperature dependence of the magnetization and susceptibility 
is identical in the two cases for reasons that are explained in the text. 




Temperature [K] 

FIG. 8: (color online) Magnetization curves and sublattice suscepti- 
bility in the traditional and extended model. The notation is as fol- 
lows: (M) and {/i} stand for the temperature-dependent magnetiza- 
tion of the Mn and Ni sublattices, respectively; stands for the 
thermal average of the absolute value of the Ni local moment; XMn.Ni 
for the sublattice susceptibilities of Mn and Ni; and A/^li„ f^; for the 
ground-state local moments, (a) Traditional model (only transverse 
fluctuations of the moment) with the interactions derived by force- 
theorem from spin-spiral calculations. The Mn magnetization shows 
a typical ferromagnetic behavior with a susceptibility peak at Tc. 
On the other hand the Ni magnetization drops rather fast, with an 
atypical susceptibility showing a plateau over a wide region, (b) 
Extended model (longitudinal fluctuations of the Ni moment are al- 
lowed) with bare Mn-Ni and Mn-Mn interactions. The Mn magne- 
tization hardly changes compared to (a), and the Tc is very simi- 
lar, but the behavior of the Ni magnetization and susceptibility are 
completely different, following the critical behavior of the Mn sub- 
lattice. At r = 0, the Ni susceptibility does not vanish, behaving 
as it is expected for a system with non-rigid moments, reaching the 
value X = l/(2a) = 0.009 /.i|/mRyd at T = 0. The average 
of the magnitude of the local Ni moment increases with tempera- 
ture. The Monte Carlo simulations included 2048 unit cells (4096 
atoms). The Ni-related quantities have been scaled up by the anno- 
tated factors in order to demonstrate the agreement with Eqs. l l38t 
and imi: (fj.) and {\n\) by a factor M^^/M^i = 14.2 and XNi by 
(MSin/Mg;)^ = 201. In particular the blue squares show the valid- 
ity of Eq. ( 142b . i.e., the derivation of the Ni from the Mn susceptibil- 
ity. The latter coincides with the XMn of the traditional renormalized 
model, as is demonstrated in Fig. [7] (bottom panel). 



Next we present a series of Monte-Carlo-simulation results. 



examining the effect and importance of the bare interactions. 
We performed simulations in the framework of the traditional 
model (only transverse fluctuations allowed) as well as the 
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extended model (longitudinal fluctuations also allowed), with 
the renormalized and bare parameters. In all cases, the Ni- 
Ni interaction was neglected, the Mn-Ni interaction was ne- 
glected for distances larger than the nearest neighbors, and the 
Mn-Mn interaction was included up to a distance of three lat- 
tice parameters. The first set of results is contained in Fig. |2l 
where magnetization curves are shown, calculated within dif- 
ferent assumptions. Here we omit showing the susceptibiliy, 
but the Curie temperature can be recognized by the charac- 
teristic inflection point of the magnetization curve. The cen- 
tral results are included in curves (a) and (e) of the top panel, 
as well as in the bottom panel, while curves (b), (c) and (d) 
merely show that neglecting some of the degrees of freedom 
or some of the bare constants leads to almost arbitrary re- 
sults. In the simulations including longitudinal moment fluc- 
tuations, more Monte Carlo steps were necessary, typically by 
an order of magnitude, in order to arrive at the same quality 
as in the simulations including only transverse fluctuations. 

Curve (a) shows the result of the traditional model with 
spin-spiral-derived (i.e., renormalized) Mn-Mn exchange pa- 
rameters, including spin-spiral-derived Mn-Ni interactions. 
This results in Tq = 870 K. Curve (b) includes the bare 
Mn-Ni parameters, but stays within the traditional model. Tq 
increases by a modest amount of 60 K, since the bare Mn-Ni 
parameters are stronger than the spin-spiral-calculated ones. 
Curve (c) also stays within the traditional model, but includes 
the bare Mn-Mn interaction, which is weaker by a factor 1/3 
compared to the renormalized value; Tc drops significantly 
to 730 K. This coincides with the experimental value, but the 
coincidence is probably fortuitous: the longitudinal fluctua- 
tions at Ni, that are essential to the derivation of the bare pa- 
rameters, are yet unaccounted for in the simulation (see, how- 
ever, the discussion on the phase space measure in subsection 
rvb . Next comes curve (d), where the longitudinal fluctuations 
are allowed for in the simulation, but taking the bare Mn-Ni 
and the renormalized Mn-Mn exchange. The increase in Tq 
(1090 K) with respect to curve (b) (same parameters but rigid 
Ni moment) is striking. The difference stems from the fact that 
the Ni local moment becomes larger at high temperatures (see 
discussion on Fig.[8]below), thus the ferromagnetic contribu- 
tion Jun-NiM ■ fjL to the Hamiltonian is effectively strength- 
ened. Finally, if we account also for the (weaker) bare Mn-Mn 
coupling in the extended model, we obtain curve (e), which is 
very close to the original curve (a), also with a very similar 
Tq ~ 850 K. This value, however, is the same that one ob- 
tains if the traditional model is used, but with the Ni moments 
completely neglected. In fact, the corresponding magnetiza- 
tion curve falls exactly on top of the Mn-contribution to curve 
(e), and the same is true for the Mn-sublattice susceptibility. 
This striking agreement is also demonstrated in the bottom 
panel of Fig. |2l and is expected on the basis of Eq. dJTl i and 
the discussion thereafter. 

Let us consider now the contrast between the traditional 
model with spin-spiral calculated interactions and the ex- 
tended model with bare interactions. Monte-Carlo results on 
these are shown in Fig. [8^ and b, representing the traditional 
and extended model respectively. The main difference lies in 
the behavior of the Ni-sublattice magnetization (/i^) and sus- 



ceptibility XNi- For better comparison we have scaled up these 
quantities. In Fig. [8^ we see that the traditional model results 
in a Ni magnetization that drops rather fast with temperature, 
much faster than the Mn magnetization. This is due to the 
relatively weak coupling between the Mn and Ni sublattices. 
Although the difference in energy scale is not directly obvi- 
ous from Fig.[3j recall that a Mn atom is surrounded by 4 only 
Ni atoms but by 12 Mn atoms at distance ai^t/V^ and 6 Mn 
atoms at distance aiat, where the exchange coupling is still 
appreciable. The Ni susceptibility, XNi shows a large plateau 
below Tq and follows the critical peak of Mn at Tc rather 
weakly. The behavior of the Mn magnetization and suscepti- 
bility, on the other hand, is characteristic of a ferromagnetic 
phase transition. 

In Fig.[8}5 we see results within the extended model. Here, 
the Ni magnetization follows the behavior of Mn. The sus- 
ceptibility starts off at a finite value at T = 0, which 
coincides with the value given by Eq. [17] in contrast to the 
vanishing susceptibility at T = for rigid-moment systems. 
At Tc the Ni susceptibility shows a peak following the crit- 
ical behavior. The Ni magnetization and susceptibility have 
been scaled up by appropriate factors to show that Eqs.(l38b 
and (l42l) are reproduced by the Monte Carlo simulation. A 
technical point to be mentioned for accuracy is that, in cal- 
culating (/i-) and {fiD in the simulation one should take the 
projection of J^i l^i the direction of the total moment at 
each Monte Carlo step, instead of making the approximation 
to calculate the average and variance of | /x; |. 

It is also interesting to see that the average magnitude of 
the local Ni moment (dotted line) increases monotonically 
with temperature. This effect, pointed out by Sandratskii^ in 
an analytical low-temperature approximation for NiMnSb, is 
connected to the fact that even above Tc there is some short- 
range order in the system, so that the equilibrium value of 
(/x) = ^Nc{M) is non- vanishing. On top of this, the fluctu- 
ations of fi shift the value of = ((/i^ + fiy + ^^)^/^) 
to higher values. The increase is expected to cease when 
fi reaches such high values that the approximation 6 = 
is no more valid (so that the fluctuations are moderated by 
the fourth-order term); this correction does not apply for 
NiMnSb, however, at least at temperatures as high as Tc, 
since the constrained-DFT calculation shown in Fig. |6] (top 
panel) yields a quadratic energy dependence at values of /i 
comparable to the ones close to Tc in Fig.[8j3. 

Experimentally, the temperature dependent magnetization 
of NiMnSb foflows Bloch's T^/^ law up to about 70- 
100 Kr^i^ which cannot be reproduced within classical 
Heisenberg models. Element (or sublattice) specific ex- 
periments were done using neutron scattering^ and X-ray 
magnetic circular dichroism (XMCD)fJ2 Neutron scattering^ 
at 15 K and 260 K shows a thermally stable Ni moment 
(0.18 /is), but a decreasing Mn moment (dropping from 
3.79 /i_B at 15 K to 3.55 /is at 260 K). Contrary to this, ac- 
cording to the XMCD results of Ref. iH] both the Ni and 
Mn moments drop rapidly at 80 K to half their ground-state 
values, and then level off up to at least 250 K. Such behav- 
ior would be counter-intuitive; the authors in |:50] write that 
surface effects possibly complicate the interpretation of MCD 
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data. 



V. FINAL REMARKS 

A. Remarks on the treatment of the weak moments and on the 
concept of renormaUzation 

Recently, Wysocki, Glasbrenner and Belashchenko^ 
(WGB) presented a study of a classical spin-fluctuation 
model. Their model Hamiltonian is similar to the one that 
we use here, with the difference that in the WGB paper all 
atoms can change the magnitude of their moments, the 4th- 
order term is not neglected, and in practice only one atomic 
species is considered in their calculations. WGB point out 
that the magnetic moments are not canonical variables, there- 
fore there is no obvious way how to choose the phase space 
measure, which should therefore be given as part of the model 
together with the Hamiltonian. In the present work we have 
chosen what they call uniform phase space measure, which 
basically amounts to dividing the /x-space in equal-volume 
infinitesimal cells with equal integration weight, and which 
is the most common choice in the literaturci^ It also amounts 
to taking a simple d^jj. integration with no further weight in 
Eqs. ( I32I36I43T |. Different choices of measure can lead to 
qualitatively different results, e.g. a fast drop of magneti- 
zation in the weakly-magnetic sublattice, as also shown by 
Sandratskii,— or even a first-order transition, as WGB find.— 
The correct measure can only be certified by the best clas- 
sical approximation to the full quantum-mechanical solution, 
which, however, remains an open problem. However, if it is 
assumed that the correct measure is not uniform but e.g. pro- 
portional to /i^^ (this was one alternative choise by WGB), 
then the weak-sublattice magnetization will drop fast at low 
temperatures, so that close to Tq only the bare parameters of 
the strong sublattice would be relevant. 

Mryasov et alM have discussed the idea of renormaUzed 
interactions in the case of FePt alloys. Our idea of renormal- 
ization is along the same lines even if not identical. The main 
novelty of the present study in this respect is the proof that the 
thermodynamic quantities can also be derived by using the 
renormalized interactions, under the assumption of a uniform 
phase space measure and of a quadratic on-site term. The lat- 
ter seems to hold true e.g. for NiMnSb and for FePt,i^2i22i^ but 
not, for example, in the case of FeRh^^il^ where higher-order 
corrections are necessary. 

Polesya et ali^ adopt a model for FePd and CoPt where 
the weak moments of Pd or Pt are determined from the strong 
moments of their neighbors via the susceptibility, similarly as 
we discuss here. However, even at higher temperatures the 
weak moments are not treated as independent variables that 
can fluctuate independently (either in direction or in length) 
but rather as enslaved quantities to their immediate neighbour- 
hood; i.e., their role in the thermodynamics is only to mediate 
an additional interaction between the strong moments. 

Bruno^"* also introduces a concept of renormalized ex- 
change parameters. However, his approach encapsulates dif- 
ferent physics than our present approach. Bruno's renormal- 



ization corrects for a systematic error, mainly due to the dif- 
ference of the constraining-field direction to the resulting mo- 
ment direction when the force theorem is applied. Our renor- 
malization on the other hand concerns the error due to the re- 
duction of the weak-moment magnitude when the strong mo- 
ments are tilted. 

Yet another concept of renormalization is described by Lou- 
nis and Dederichs.^- Using a multiple-scattering approach, 
they consider the energy expansion as a function of the angles 
between moments. As they find, at high angles corrections 
are necessary to the phenomenological Hamiltonian (e.g., bi- 
quadratic or four-spin terms). However, at low angles these 
corrections can be partly included in the Heisenberg model 
via a renormalization of the exchange parameters, recovering 
correct energy scales and Curie temperatures. 

B. Remarks on the prediction of the Curie temperature 

The Curie temperatures calculated within the mainstream 
approach to the adiabatic spin-dynamics are in many cases in 
agreement with experiment to within 10-15%, but with no ob- 
vious systematics toward over- or underestimation. The main 
source of error is not clear Considering the most serious ap- 
proximations made, error can stem from 

• (i) the use of local density functional theory (LSDA or 
GGA) for total energy calculations without further cor- 
rections for electron exchange and correlation. 

• (ii) the use of the adiabatic approximation, 

• (iii) the assumption of a classical, rather than quantum, 
Heisenberg model, 

• (iv) the assumption that the exchange constants do not 
change as a function of temperature, 

• (v) the assumption of rigid spins of the strong moments 
in the Heisenberg model. 

In general, these factors have possibly different weight for dif- 
ferent materials. Concerning point (i), theories that provide a 
better treatment of correlations exist, e.g., the LSDA-nt/ or 
LSDA combined with dynamical mean field theory at zero or 
finite temperatures. Also within such theories exchange pa- 
rameters can be derived (see, e.g., Ref.flTl). However, param- 
eters are required (as is the Coulomb repulsion U), and it is 
usually not obvious how to determine these uniquely. 

As for point (ii), there are promising developments in the 
calculation of magnon spectra (including magnon lifetime ef- 
fects) within time-dependent density functional theory^^— or 
many-body perturbation theory^. These can prove very use- 
ful in the future (they can also be extended to finite tempera- 
tures), however at this point they are computationally too de- 
manding for systematic calculations of the Curie temperature. 

Point (iii) (the classical assumption), can be improved upon 
by using the random phase approximation to solve the quan- 
tum Heisenberg model. However, for itinerant electron sys- 
tems the local moment M does not corespond to some in- 
teger or half-integer value of the spin S either in the form 
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M = y^S{S + 1) or in the form = S' ^b- In fact, cal- 
culations of Heusler compounds in Ref. 0] have shown that 
for a reasonable choice of S the Curie temperature is strongly 
overestimated, while the classical limit of the random phase 
approximation, with a choice of large S (with appropriate 
normalization of the exchange parameters so that the prod- 
uct JijSiSj remains constant), results in reasonable values of 
Tq- Therefore, a quantum Heisenberg model is perhaps a bet- 
ter for a correct description of the shape of the magnetization 
curve Af(T), but a poor choice for a correct Tc, at least in 
itinerant electron systems, if the exchange constants are cal- 
culated within the adiabatic approximation. This conclusion 
is in accordance to the spirit of adiabatic spin dynamcics,'- 
where the effective interactions ,/„„/ correspond to the equa- 
tion of motion of the expectation value of the local moments, 
i.e., to classical quantities, not operators. 

Corrections to point (iv) can be treated within local den- 
sity functional theory if the exchange constants are calculated 
starting from a disordered local moment state. This requires 
use of the coherent-potential approximation (CPA), and has 
been proposed for example in Ref. |3]. The use of the CPA 
for the description of the disordered local moment state at 
Tc underestimates the existence of magnetic short range or- 
der (which is known to be present); however, it constitutes a 
promissing approach, since it can be systematically improved 
e.g. by the use of a non-local CPA.^ 

Finally, point (v) becomes a serious approximation in sys- 
tems of weak magnetic moments, such as ferromagnetic Ni, 
and has been widely discussed in the literature as we noted 
in SecHVl Corresponding corrections for multiple-scattering 
based methods have been recently proposed, e.g., by Bruno^ 
and Shallcross and co-workerSf^ In a more recent work by 
Ruban et al.r^ based on an expansion of the energy within 
the disordered local moment state, promising results were ob- 
tained showing the fundamental importance of longitudinal 
corrections to the local moment for Tq in ferromagnetic Ni. 



VI. SUMMARY AND CONCLUSIONS 

In the first part of this work we have presented a way to cal- 
culate the interatomic exchange constants within the FLAPW 
method for multi-atomic lattices. It is based on the concept 
of adiabatic spin dynamics. The exchange constants are har- 
vested by an inverse Fourier transformation involving static 
spin-spiral energies. Symmetry relations obeyed by the spin 
spiral energies have been found to greatly reduce the numer- 
ical effort, in particular regarding confinement of the inverse 
Fourier transformation in the irreducible wedge of the Bril- 
louin zone. Furthermore, the force-theorem approximation 
has been tested and found to be adequate for small cone an- 
gles of the spin spirals. However, we have shown that appli- 
cation of the force theorem requires special treatment of the 
intersitial region, namely setting there the magnetic part of 
the exchange-correlation field to zero. 

In the second part of the present work we have proposed 
a way to explore multicomponent systems where a magneti- 
cally strong sublattice coexists with a magnetically weak sub- 



lattice, necessiating a consideration of longitudinal changes 
of the local magnetic moment. We find that, under the fre- 
quently met condition of a parabolic dependence of the en- 
ergy on the weak-moment magnitude, the weak moments and 
their interactions can be eliminated via an analytical integra- 
tion of the partition function in favour of the strong moments 
with renormalized, temperature-independent exchange con- 
stants. We also show that the renormalized constants are actu- 
ally the ones probed by constrained spin-spiral calculations 
of the strong-moment subsystem, simplifying calculations. 
This renormalization will affect various quantities such as 
temperature-dependent magnetization, susceptibility, or spin- 
stiffness constant. 

Our calculations of the exchange constants and the Curie 
temperature of fee Co, NiMnSb, and Co2MnSi have shown 
reasonable agreement with experiment, with the calculated 
values for Tc being within 15% of the experimentally ob- 
tained results. However, there is no systematic trend in the 
theoretical error; the experimentally obtained values can be 
underestimated in some cases (e.g., fee Co), overestimated in 
others (e.g., NiMnSb), or even reproduced with good accuracy 
(e.g., CozMnSi). 
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Appendix 

We provide a derivation of a formula concerning the renor- 
malization of interactions J^j of the strong-moment subsys- 
tem, {Mi},i G {1, • • • ,Ns} when it is in contact with a 
weak-moment subsystem, {fii}, I, I' G {1, • • • , iV^}, in the 
presence of interactions Aw among the moments of the weak 
subsystem. The main complication in the presence of A^' ^ 
for I ^ I' is that the strong moments are interacting with a sys- 
tem of coupled harmonic terms, instead of independent har- 
monic terms which were treated in Sec. |IV] In particular we 
adopt the following conventions. The Hamiltonian reads 



n = ns + n„ + Hint (44) 
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where the strong-system, weak-system and interacting parts 
are respectively: 



= -Wjt,M,-M, 



Hint 



(45) 

(46) 
(47) 



In Eq. (l46b , the diagonal part An = a > is the quadratic 
on-site energy term while the off-diagonal terms Aw describe 
intersite interactions between the weak momentsi^ In case of 
ferromagnetic coupling it is expected that Aw < for I ^ 
(but with the determinant det| A;;/ 1 > 0). 



J 



The strategy is to eliminate the weak moments by inte- 
grating analytically the weak plus intercting part of the parti- 
tion function, ending up with renormalized interactions of the 
strong moments. To this end we take advantage of the identity 



dxi ■ ■ ■ dxN cxp 



E 



XX' 



rN/2 



exp 



Y.iC-')xx'bxbx' 



XX' 



(48) 



where det C is the determinant of the positive-definite ma- 
trix C. It is convenient to use a combined index A ~ {l,a) 



with a e {xyz} and define b\ 



1 



-T 



n{l) 



Mia and 



Cxx' = j^AwSac.'- Then iC-')xx' = kBT{A-^)w6a 
Applying this to the partition function Z yields 



(detA)3/2 
(detA)3/2 



dfli ■ ■ ■ d^N^ exp 
drii ■ ■ ■ drijVa exp 



1 



1 



l~Lr 



cxp 



1 / K \2 



knT 



© E(^-^)«'EE^"-^"' 



n{l) n(l') 



where a renormalized Hamiltonian of the strong-moment sublattice has been introduced, 

^- = E • + (^)' E(^"')«' E E • 



n{l) n'(l') 



(49) 
(50) 

(51) 
(52) 



This reduces to Hr of Eq. OTI ) if Aw — adw- Expres- 
sion ( l52b is obviously of the traditional Heisenberg type, but 
a further reduction to a form with renormalized parameters, 
Hr = — ^ J'ij^i ■ requires knowledge of the specific 
geometry of each problem taking into account the sums over 
neighbours of I and V , J2n{i) ^"^^ J2n'{i')- 

Thermal averages (/({TVfi})) can be calculated by 



(/({M.})) 



1 (TTkBT) 



3JVw/2 



dr!i---dl)jvJ({M,})cxp 



(det A)3/2 



knT 



(53) 



thus the factor 



cancels and the determinant dot A 



(dctyl)3/2 

need not be calculated. Thus one ends up with a usual 
Heisenberg-model treatment. To gain some more insight, one 
can recognize that in the presence of non-diagonal Aw we 
have, formally, a system of coupled harmonic oscillators in- 
teracting with the strong moments. The normal modes of the 
coupled oscillators are itinerant, and therefore the renormal- 
ized interactions are also long-ranged. Even if Aw are short 



ranged, making the matrix A sparse, Eq. ( |52] | shows that the 
renormalized interactions involve the matrix A^^, which is 
normally not sparse. 



To complete the circle, it has to be shown for practical 
applications that the renormalized interactions appearing in 
Eq. ( |52] | are the quantities that are probed by a DFT calcu- 
lation where the directions of the strong moments are con- 
strained. In other words, in such a DFT calculation the weak 
moments are allowed to relax to their equilibrium values un- 
der the directional constraint on the strong moments. The 
question is, if then the energy dependence ( l52l i is recov- 
ered, assuming that Eqs. ( |45p7l ) are a good approximation 
to the DFT energy landscape. The answer is straightforward 
if we calculate the total energy of the constrained Heisenberg 
Hamiltonian, i.e., without an integration over the Af j. The 
constrained partition function is just the last term of Eq. ( |49l ), 
Zi ^ /d^^i • • •(i3^^^e-('""+^'"')/*=«^. The total energy 
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at T ^ is 



H.({M,}) 
H.({M,}) 



kT^ dZi 



Zi dT 



(54) 



T=0 



2^) E(^"^VEE^"-^"' (55) 

n{l)n'{l') 



q.e.d. This means that constrained (spin spiral) DFT cal- 
culations on the strong sublattice are already corresponding 
to the renormalized Hamiltonian and are therefore by virtue 
of Eqs. ( 14915 II ) sufficient for the calculation of the strong- 
sublattice thermal averages, without the need to calculate the 
bare parameters or An . 

However, if the weak-sublattice thermal averages are to be 
calculated, one must additionally gain knowledge on the ma- 
trix Aw of Eq. ( |46] l as well as of the bare parameters 
of Eq. ( |46] | and of k in Eq. (l47b . The scheme presented in 
Sec. IIV Al shows how this can be done in the rather simple 
case of NiMnSb (where An = aSn is diagonal), but in gen- 
eral this problem must be solved according to the geometry 
and other factors in each case. In particular for the calcula- 
tion of the off-diagonal Aw , probably it is easiest to calculate 
directly the susceptibility matrix {A~^)w by applying in the 
DFT calculation a longitudinal external field on one atom and 



probing the response in the moment of the neighboring atoms, 
or to make a transformation to the normal modes in Fourier 
space and calculate A{q). 

If all the ingredients are available, then one can calculate 
thermal averages and correlation functions of the weak mo- 
ments either by a direct Monte-Carlo calculation, or by reduc- 
ing these correlation functions to correlations of the strong 
moments in the presence of only the renormalized Hamilto- 
nian, in an analogous way to the case of Eq. ( l43T i. where the 
parameters Aw and k must be inevitably contained in the 
expansion coefficients. We present an outline of how this 
is achieved in practice. The complication here compared to 
Eq. (03]) is that the matrix Aw is not diagonal. It is, how- 
ever, real and symmetric, therefore one proceeds by bringing 
it to a diagonal form D with elements Di. If R is the di- 
agonalization matrix, then D = RDR^, and we define the 
transformation of the moment-coordinates figa ~ J2i Rqil^ia 
and Mqa = '}2,,^RqnMna- Also sincc Aw is real and sym- 
metric the phase-space element is unchanged: (Pjl = dPfi. 
(This transformation is analogous to one bringing a system 
of coupled harmonic oscillators in a normal-mode represen- 
tation.) Then one has, for the arbitrary correlation func- 
tion among the xyz-components of the A^w weak moments, 
((Mia;)™'" • • • (/"Afwz)™"~">, the following integral: 



J 



la k I 



-Diiili 



-^AwHla 



U'a 



J d^iii ■ ■ ■ d'fiNjY. RikP^-T''' • • • (E Rl^kf^^^^T"-' exp J ^ 



2a 



EE^^' 

loL n{l) 



ksT 

1 K 

A 2^ 



-EamL + I-EE^^' 



la 



2a 

let n(l) 
1 , K 



naf^la 



2(1) 



(56) 



(57) 



(58) 



The last expression can be handled by expanding the term 
(^^ Rf^.p.ka)^'''" and proceeding to an analytic integration of 
powers of /i times the exponential just as in Eq. i43[ . Mixed 



correlation functions between strong and weak moments can 
also be reduced to strong-moment correlation functions in the 
renormalized model via this prescribed route. 
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